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Abstract 



In this article we have successfully obtained an exact solution of a particular case 
of SIR and SIS epidemic models given by Kermack and Mckendrick Jl^ for constant 
population, which are described by coupled nonlinear differential equations. Our re- 
sult has no limiting conditions for any parameter involved in the given models. In 
epidemiology many researchers believe that it is very hard to get an exact solution for 
such models. We hope this solution will be an opening window and good addition in 
the area of epidemiology. 
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1 Introduction 

The mathematical modelling of spread of diseases is an area of biology named as epi- 
demiology. Physical real life problems arising in epidemiology may be described, in a 
first formulation, using differential equations. As far as our knowledge is concerned 
the first mathematical model of epidemiology was formulated and solved by Daniel 
Bernoulli in 1760. Since the time of Kermack and McKendrick [T], the study of math- 
ematical epidemiology has grown rapidly, with a large variety of models having been 
formulated and applied to infectious diseases [2]-|l]- 



The dilemma with models in epidemiology is, sometimes even the simplest mathemat- 
ical model of natural phenomena with sets of first-order ordinary differential equations 
are non-integrable. Nucci and Leach [5] have discussed these features and obtained 
an integrable SIS model by applying Lie analysis. According to Nucci and Leach [5] 
It would seem that a fatal disease which this models is also not good for mathemat- 
ics!. Khan et.al., [6] has obtained the convergent series solutions for SIR and SIS 
models using homotopy analysis method. Considerable attention has been given to 
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these models by several authors using stability, bifurcation theory, Lyapunov method, 
Poincare-Bendixson type theorems, index and topological concepts see ([2]-[T6]). 

To understand what SIS and SIR stands for let us consider some constant population, 
let divide them into three components: the susceptible component, S, who can catch 
the disease; the infective component, I, who are infected and can transmit the disease 
to susceptible, and the removed component, R, who either had the disease and recov- 
ered, died or have developed immunity, or have been removed from contact with other 
components. SIR model contains all three classes, while, in SIS model the infectives 
return to the susceptible class on recovery because the disease confers no immunity 
against reinfection. Such models are more effective for disease caused by bacteria or 
helminth agents and also for most sexually transmitted disease. 

In fact closed-form solutions to mathematical models play a vital role in the proper 
understanding of qualitative features of natural science. In epidemiology mathemati- 
cal models describes some effects of disease in some population so one need to be very 
careful while analyzing such models. Keeping this view of point, here in this article 
our main goal is to obtain pure analytic exact solution to the SIS and SIR epidemic 
model which is a particular case of the model given by Kermack and McKendrick [1] . 

In second section we present the solution for SIS model while in third section we 
present solution for SIR model. 



2 Mathematical Analysis for SIS model 

The SIS model by Kermack and Mckendrick [T] is given as follows 

s'{t) = —rsi + ai, (1) 
i'{t) = rsi — ai, (2) 

subject to Initial conditions, 

z(0) = /o, s(0) = 5o, (3) 

where r > 0, Jq > and 5*0 > 0. Here r is the infectivity coefficient of the typical 
Lotka-Volterra interaction term and a the recovery coefficient. From equation ([1]) 
and (|2]) one can easily observe that s{t) + i{t) = k, where k is total population. This 
model is different from SIR model as in this model the recovered members return to 
the class S at a rate of ia instead of moving to class R. 
Now using s{t) = k — i{t) in equation ([2]), 



i' = r {k — i)i — ai. 
Integrating the above equation (jl]) with substitution y = i^^, we get 



(4) 
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{rk-a)t ^ / ^^irk-a)t^^ ^ 



ye' 



where C is integrating constant. After simplifying equation. ([5]) and letting (3 
rk — a, we get, 



y = ^ ■ (6) 



Using back substitution y = i ^, 



r + ^Ce- 



After applying initial condition z(0) = iq on equation ([7]) we get, 

^ ^ /3 - ipr 

/3zo 

The exact solution for s{t) obtained by using the fact s{t) + i(t) = k is, 

/3 



s{t) = k- 



3 Mathematical Analysis for SIR model 

The SIR model by Kermack and Mckendrick [IJ is given as follows 

s'{t) = -l3s{t)i{t) - fis{t) + ij, (8) 
i'{t) = I3s{t)i{t) - iii{t), (9) 

subject to Initial conditions, 

z(0)=zo, s(0) = so, (10) 

(') denotes the derivative with respect to time, where /3 > 0, io > and Sq > 0. Here 
13 is the infectivity coefficient of the typical Lotka-Volterra interaction term and the 
recovery coefficient. Now adding equation ([8]) and equation (j9]). 



(s + i)' = -fi {s + i) + /i. 
Integrating the above equation [11], we get 

s{t) = l + Ce(-^)*-z(t), 



(11) 



(12) 
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where C is an integrating constant. Now using equation (fT2|) in equation ([9]), we get, 

i'{t) = {l3-fi + f3Ce'^'')i{t) - (3i\t). (13) 
Integrating the above equation with substitution z = we get 



Where D is an integrating constant. Now observing that e"'^* approaches to zero as 
t ^ oo, here we take series expansion for e~'^* and neglecting square and higher terms 
we get e"'^* = 1 — fit. Substituting back in equation (fT^ one gets, 



^g(/3-M)t-(^)(l-Mt) 



Z 



J Pe^('-^^'-^'^^^'-^'^dt + D, (15) 

1^ . T^.-(8-u+BC)t^^ 



^g-(/3-M+/3C)igV, (16) 



(3-fi + (3C 

Substituting z = i^^ and /3 — /i + /3C = A, in equation (fT6!) we get 

= ^ jc- (17) 

Applying initial condition z(0) = io and s(0) = Sq in equation (fT2!) and (!20!) we get, 

A - 2o/3 



D 



C = So + to-l. (18) 
The exact solution for s{t) obtained by using the value of i{t) in equation (fT2l) which 



IS, 



s{t) = 1 + (so + - 1)(1 - ^t) -— ^ (19) 



Xige M 



'(t) = ,,, .(.o+.n-i) - (20) 

/3 + A(^gfer)e-^*+^^ 

Aioe A* 



4 Conclusion 



In this paper an exact analytic solution to a particular case of SIR and SIS models 
with constant population is obtained with the help of direct integration tools. As these 
models have already been observed by different researchers with different techniques 
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but we know the importance of exact solution, that exact solutions allow researchers 
to run experiments, by introducing natural initial or boundary conditions, to observe 
the behavior of functions involved in equations. This present work will be of great sup- 
port to the epidemiology to investigate the hidden phenomenon and validity of these 
models which is already being considered and used by different researchers. This kind 
of approach will attract other nonlinear models as well having been analyzed again 
for the possible existence of exact solution. 
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